On the possible role of massive neutrinos in cosmological structure formation 



Massimiliano Lattanzi, 1, 2 Remo Ruffini, 1 ' 2 and Gregory Vereshchagin 1,3 

1 ICRA — International Center for Relativistic Astrophysics. 
2 Dipartimento di Fisica, Universita di Roma "La Sapienza" , Piazzale Aldo Moro 5, 1-00185 Roma, Italy. 
3 Byelorussian State University, Theoretical Physics Department, 
Skorina ave. 4> 220050, Minsk, Republic of Belarus. 

In addition to the problem of galaxy formation, one of the greatest open questions of cosmology 
is represented by the existence of an asymmetry between matter and antimatter in the baryonic 
component of the Universe. We believe that a net lepton number for the three neutrino species can 
be used to understand this asymmetry. This also implies an asymmetry in the matter-antimatter 
component of the leptons. The existence of a nonnull lepton number for the neutrinos can easily 
explain a cosmological abundance of neutrinos consistent with the one needed to explain both the 
rotation curves of galaxies and the flatness of the Universe. Some propedeutic results are presented 
in order to attack this problem. 
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I. EVIDENCE FOR DARK MATTER AND THE POSSIBLE ROLE OF NEUTRINOS 

The most popular model of the Universe being currently discussed in the literature is one usually indicated by the 
"ACDM model" (see e.g. E.W. Kolb, in these same proceedings [lj), which implies that almost 75% of the energy 
density of the Universe is due to a cosmological term, while only 1% is due to neutrinos This result seems 

to be in conflict with our current knowledge of the rest of physics. We will give in the following some propedeutical 
considerations to reconsider this problem. 

Many theoretical considerations and observational facts make it clear (see and references therein) that luminous 
matter alone cannot account for the whole matter content of the Universe. Among them there are the considerations 
on cosmological nucleosynthesis as well as the measurements of the the cosmic background radiation (CBR) 

anisotropy spectrum ^?|~[23- In both cases the fit is consistent with a cosmological model in which just a fraction 
smaller than 10% of the total density is due to baryons. 

Strong evidence for the presence of dark matter is directly given by the rotation curves of galaxies pH ] . If we assume 
for simplicity a spherical or ellipsoidal mass distribution for a galaxy, the orbital velocity at a radius r is given by 
Newton's equation of motion: 

J = «, (1) 
r 

where G is the gravitational constant and M(r) is the mass contained in a sphere of radius r. The peculiar velocity of 
stars beyond the visible edge of the galaxy should then decrease as 1/r. What is observed instead is that the velocity 
stays nearly constant with increasing r. This requires a halo of invisible dark matter to be present outside the radius 
of the visible matter. From observations it follows that the halo radius can be 10 times larger than the radius of 
visible part of the galaxy. Then from Eq. it follows that M^aio is at least 10 times larger than the galactic mass 
Mg a i. We now assume that galactic halos are composed of neutral fermions of mass m x and apply a gravitational 
Thomas-Fermi model to this system. The equation for the dimensionless gravitational potential ^ as a function of a 
radial coordinate x is 

d 2 X X 3/2 

dx 2 yfx ' ^ 

with the initial condition x(0) = 0. For any number of particles N this equation can be used to compute the radius 
R of the system. Since the total mass is M = Nm x , this defines a relation between M, R and m x which allows 
one to estimate m x , since both the value of the total mass and radius are known. Using for the radius of the halo 
Rhaio = 10-Rgai and for the mass of the halo Mhaio — 10M SO ; we obtain the particle mass m x « 4eV. This value, 
although obtained by a simple argument, is in good agreement with the recent constraints on the electron neutrino 
mass obtained from the spectrum of tritium beta decay (m„ < 2.5 eV H^)- 

This agreement can even be improved if one takes into account different families of neutrinos. In this case, the 

1 /4 

effect of the additional degrees of freedom can be expressed in terms of the effective mass m e // = (X^ILi m f) ■ If 
we consider three families with nearly the same mass, they provide a factor 3 1 / 4 ~ 1.32 to improve the bound to 
m„ ~ 3 eV. 

Neutrinos were considered as the best candidate for dark matter about thirty years ago. Indeed, it was shown by 
Gerstein and Zel'dovich in 1966 23] that if these particles have a small mass m v ~ 30 eV, they provide a large energy 
density contribution up to the value Q. v ~ 1. It is in effect easy to show that the density parameter of a single family 
of mass m„ <C 1 Mev would be, assuming a null chemical potential: 

°* = < 3 > 

where h is the Hubble constant in units of 100 km sec -1 Mpc -1 . The generalization to the case of a nonnull chemical 
potential is (see e.g. Rumni & Song |24j): 



(4) 



where £ is the dimensionless chemical potential (so called degeneracy parameter) at decoupling and 
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|e| 3 +4r ? (2)|e|+4^(-l) 



fc+i £ 



(•5) 



Here 77(71) denotes the Riemann r\ function of the index n. This allows to obtain even higher values of Q„. In fact we 
have carried out a more detailed analysis, taking into account three different neutrino flavours with the same mass but 
different chemical potentials. Using the limits that primordial nucleosynthesis imposes on the degeneracy parameters 
of the e, [i and r neutrinos ^(|, we found that if m„ ~ 2 eV and £ e ~ 0.4 then VL U ~ 1 |25j . 

However, in 1979 Tremaine and Gunn prj claimed that massive neutrinos cannot be considered as a dark matter 
candidate. Their paper was very influential and turned most cosmologists away from neutrinos as cosmologically 
important particles. In their paper, Tremaine and Gunn establish an upper limit to the neutrino mass of m„ < 1.2 eV. 
They obtain this limit from arguments based on the velocity dispersion within galaxies and clusters, under the 
hypothesis that m„ < 1 MeV, so that neutrinos are ultrarelativistic at decoupling. They also obtain a lower bound 
m v > 20 eV from an argument based on phase space density considerations and on the rotation curves of galaxies. 
In this very strange situation with wildly contradictory constraints, they possibly see a way out by avoiding the fact 
that neutrinos are ultrarelativistic at decoupling and they conclude that massive galactic halos cannot be composed 
of stable neutral leptons of mass < 1 MeV. 

While the Gunn and Tremaine result deeply influenced astrophysicists against the possible role neutrinos in cos- 
mology, especially in the U.S.A., in 1977 Lee and Weinberg [23 turned their attention to massive neutrinos with 
m„ > 2 GeV, showing that such particles could provide a large contribution to the energy density of the Universe, 
in spite of a much smaller value of their number density. This paper was among the first to consider very massive 
particles as candidates for dark matter. This very interesting work, together with the Gunn and Tremaine purported 
difficulties for the neutrino scenario, induced some cosmologists to turn their attention to very massive particles, thus 
marking the birth of cold dark matter models. 

A clear counterexample to the Gunn and Tremaine bound, which was indeed derived from a nontransparent mixture 
of quantum limits on classical Maxwell-Boltzmann statistics, was given by Gao and Ruffini .28]. They established 
a different upper bound on the neutrino mass from the assumption that galactic halos are composed of degenerate 
neutrinos: m„ < 15 eV. This result was further developed and confirmed by Arbolino and Ruffini |29j. They explicitly 
showed that rotation curves for galaxies in perfect agreement with the observations can be obtained for neutrino masses 
of the order of 9 eV. This limit could be lowered further if semidegenerate configurations for the neutrino halo were 
to be considered (see e.g. Merafina and Ruffini [23|"|23)i Ingrosso, Merafina, Ruffini and Strafella 32] ). 

Today, quite apart from the rotation curves of galaxies, the recent determination of the neutrino masses would 
appear to be in contradiction with an assumption of Q v ~ 1. However, this is only an apparent difficulty, since for 
semidegenerate distributions (see Eq. Q) this equality can indeed be fulfilled and important consequences on the 
matter-antimatter asymmetry in the leptonic component of the Universe can be inferred. 

One of the most interesting features of neutrino cosmology is that they establish a natural cutoff for the largest 
possible structure in the Universe, related to the maximum value of the Jeans mass when the neutrinos become 
nonrelativistic: 



where z n r is the redshift at which neutrinos enter the nonrelativistic regime, g v is the number of quantum degrees of 
freedom, N„ is the number of neutrino families, A(£) is as defined in Eq. ©, and 



This mass appears to be essential in determining the upper cutoff for a possible fractal structure of the Universe. 
Real difficulties still exist today in understanding the details of the fragmentation of these masses of 10 17 M© and 
the development of smaller structures all the way down to galaxies. In this lecture we illustrate some of this basic 
problems which still need additional work before a detailed correspondence with observations can be obtained. 




(0) 




(7) 



II. LARGE SCALE STRUCTURE 



A. The cosmological principle 



There have been three distinct moments in the development of the so called cosmological principle which is at the 
very basis of our approach to the analysis of the Universe. The first formulation of the cosmological principle can be 
simply stated: 



Such a cosmological principle was enunciated a few years after the introduction of the field equations of general 
relativity by Einstein himself in the quest for visualizing a Universe the most democratic with respect to any 
special point and any possible moment of time: a Universe everlasting in time and totally homogenous in the spatial 
directions. No solution fulfilling such a cosmological principle could be found, and Einstein was so strongly confident 
of the validity of this principle that he modified his field equations of general relativity by introducing a cosmological 
constant A. George Gamow said that Einstein later on considered that the biggest mistake in his life. 

It was through the work of Alexander Alexandrovich Friedmann |34l|-|35| that a new cosmological principle was 
advanced: 



As long as we look at our 'neighbour' Universe, this statement is certainly false, because the distribution of matter is 
far from homogeneous: there are planets, stars, and going to larger scales, galaxies and clusters of galaxies separated 
by almost empty regions. However, the Friedmann principle should apply when we average this distribution over a 
volume containing a large enough number of galaxies. For such a spatially homogeneous Universe Friedmann [34[ 
found in 1922 explicit analytic solutions of the Einstein equations of general relativity. A remarkable property of 
these solutions is that they describe a non-static Universe. At that time, there was no observational evidence for the 
temporal evolution of the whole Universe. The first evidence came in 1929 from the observation by Hubble [3j| of 
the recession of the nebulae. Hubble was the first trying to study the spatial distribution of objects as large as the 
galaxies, at that time thought to be the largest self-gravitating systems to exist. The Hubble law, interpreted within 
the framework of Friedmann cosmology, implied that the galaxy distribution is close to homogenous on the large-scale 
average |37j-|40j. 

It was through the above mentioned work of Hubble together with the later remarkable work of George Gamow 
and his collaborators (1946-1949) l4ll-|H(j who postulated an initially hot Universe, and the detailed work of Fermi 
and Turkievich in the same years pl\ introducing the first computation of cosmological nucleosynthesis, that the 
Friedmann Universe has grown to become the standard paradigm in cosmology following the discovery of the CBR 
by Penzias & Wilson in 1965 |5^ . 

In effect, one of the strongest predictions of Big Bang model is the presence of a background microwave radiation, 
a relic of the early Universe. This radiation is highly isotropic, reflecting through the coupling with matter the 
high isotropy and homogeneity of the primeval plasma. This tells us that the cosmological principle, and then the 
Friedmann picture, safely apply to the early Universe. Homogeneity on very large scales is confirmed by present day 
observations of in particular: 

• X-ray background [Hflj . 



• radio sources 




• gamma ray bursts distribution |55j . 

• galaxies and clusters of galaxies . 

So much so for the very large scales, but what about structures like galaxies and clusters of galaxies? More and 
more they appear to be distributed without any apparent homogeneity, but on the contrary showing regularities in 
an ap pare nt hierarchical distribution of galaxies, clusters of galaxies and superclusters of galaxies separated by large 
voids 

[E3-IH3 (see Pig. HJ. Slowly but more and more clearly the presence of a fractal distribution in the Universe 
has started to surface and with it a new cosmological principle which can be simply expressed: 



All the events in the Universe are equivalent. 



(8) 



All the points in the Universe are equivalent. 



(9) 



Figure 1: The distribution of galaxies in the 2dFGRS (from |S^V Courtesy of J. A. Peacock and the 2dFGRS Team. 



All the observers in the Universe are equivalent. (10) 

We shall recall in the following a few basic points which are essential in arriving at this new principle and make 
possible the verification of its possible validity. 

B. Two-point correlation function 

The statistical description of clustering is based on the concept of correlation, namely, more precisely, on the 
probability of finding an object in the vicinity of another one. The standard way to quantify this probability is to 
define the two-point correlation function £(x) 

Consider a distribution of objects in space described by the number density function n(x). The probability that an 
object is found in an infinitesimal volume SV centered around the point x is proportional to the volume itself: 

SPocSV. (11) 

In the absence of structure, the joint probability of finding two objects in two different infinitesimal volumes 5Vi and 
SV2, centered respectively around x\ and X2, is given by the product of the two probabilities: 

SP = 6P 1 6P 2 oc 8Vx8V 2 . (12) 

On the other hand, if objects have a tendence to cluster, we will find an excess probability: 

SPac8Vi5V2-(l + t@i,&)). (13) 

According to the cosmological principle, we don't expect the correlation function to depend either on the position or 
on the direction, but only on separation beetween volumes: £(a?i,2?2) = £( r i2)> where r\2 = \x\ — x 2 \- 
An equivalent definition of the two-point correlation function is the following: 

£(r 12 ) =< 5{xx)5{x 2 ) >, (14) 

where < ... > denotes averaging over all pairs of points in space separated by a distance ria, and 5(0) = (n(x) — n)/n. 



C. Observed galaxy distribution 



Observational data coming from galactic surveys are usually expressed in the form of a correlation function £(tt, a) 
in redshift space, where 7r is a separation along the line of sight and a is a angular separation on the plane of the sky 
between two galaxies. It is then possible to obtain the real-space correlation function this step is never a trivial 
one, but we will not go into details since it is beyond the scope of this review. 

Peebles j6^ | has shown that the distribution of galaxies can be described by a two point correlation function with 
a simple power law form: 



, 1.77 

r 



U r ) = — - r < lOfc^Mpc, (15) 



r„ , 

where h is the present day Hubble parameter measured in units of 100 s ^ pc - The correlation length r g determines 
the typical distance between objects. For galaxies, it has been estimated to be ~ 5h~ 1 Mpc. 

For clusters of galaxies the same power law was found first by Bahcall and Soneira |r34j and then Klypin and Kopylov 

M 

£ c ( r ) = (_J , 5k- 1 < r < 150/i _1 Mpc (16) 

with different correlation lengths, namely r c ~ 25/i _1 Mpc. Furthermore, Bahcall and Burgett j6|| found a correlation 
function for superclusters of galaxies with the same power law. 

Recent observations support these conclusions. The first results from the Sloan Digital Sky Survey (SDSS) on 
galaxy clustering [67j for about 30,000 galaxies give a real-space correlation function of 

/ v -1.75±0.Q3 

£ g (r)=(I.\ , o.lh- 1 < r < 16fc -1 Mpc, (17) 

where ro ~ 6.1 ± 0.2/i _1 Mpc. The geometry of samples in SDSS is quite close to the Las Campanas Redshift Survey 
|68| and the results are very similar, but with much better resolution. 

The largest data set today is a 2dF Galaxy Redshift Survey |63 (see figQ) that consists of approximately 250,000 
galaxy redshifts. Their results are: 

/ x -1.87 

e s =f-J r ~ (6 + lO^Mpc, 

ffl=(-J r ~(3-6)^- 1 Mpc. 

Their measurements are in agreement with previous surveys. However, having much smaller statistical errors they 
were able to find a slight difference in the power law exponent as well as in the correlation length for distances or 
redshifts, colors and types of galaxies. The result can be found at [Tcj . 

D. Power law clustering and fractals 

It is clear that once a correlation function is given, the density of objects around any randomly chosen member of 
the system is: 

n(r) oc 1 +f(r). (18) 

If the correlation function has a power law behaviour with exponent 7 then: 

£(r)ocr~ 7 . (19) 

As is the case for galaxies and clusters of galaxies, where 7 ~ 1.8, the number of objects in a given volume scales in 
a similar way: 

N(r) (xr 3 -' ( . (20) 



So for noninteger 7, the number of objects scales with a fractional power of the radius of the volume under consider- 
ation. This behaviour is typical of fractal sets. 

A fractal is a set in which 'mass' and 'radius' are linked by a fractional power law |7lj : 



where Dp is the fractional or Hausdorff dimension of the set. So galaxies seem to show, at least up to scales of about 
100 Mpc, a fractal distribution with Dp ~ 1.2. 

A crucial characteristic of a fractal distribution is the presence of fluctuations at all length scales and consequently 
the impossibility of defining an average value for the density. In a simple fractal set each observer at a matter point 
belonging to the set observes the same matter distribution as any other observer belonging to the set. In this sense, 
the fractal naturally leads to the formulation of the cosmological principle expressed by the statement Ijl0|l . Three 
very serious issues, however, arise for the compatibility of a fractal structure with that of a cosmological Fricdmann 
model. (1) There must necessarily be a cutoff in the fractal distribution (see Fig. [2Jin order to recover the overrall 
homogeneity observed at the large scales up to z =~ 10 3 in the CMB. (2) Such a cutoff must occur homogeneously 
all over the Universe. (3) The dimension of that cutoff must automatically originate as the characteristic length 
determined by the microphysical properties of the dark matter composing the Universe. These were the three basic 
thoughts which motivated the cellular model of the Universe introduced by Ruffini in the eighties |z3|-|s3- ^ n this 
model fractals arise from successive fragmentation of primordial structures, the so called 'elementary cells', formed 
via gravitational instability in the neutrino component of the matter of the Universe. We shall further expand on this 
idea in the following paragraphs. 



Gravitational instability is usually considered to be the basic mechanism of structure formation in the Universe 
(see for example |15|1. It is believed that small inhomogeneities are already present at some initial time in the early 
Universe. Such small perturbations will grow due to gravitational attraction, because overdense regions will accrete 
matter from the neighbouring regions, increasing the density contrast. 

One of the simplest examples showing the process of gravitational instability is a perfect fluid model. If density 
distribution in a self-gravitating fluid is slightly nonuniform, i.e. small density perturbations exist, they will tend to 
grow. When the density contrast is small, the linear approximation can be used. The main advantage of linear theory 
is that perturbations on different scales evolve independently. 

It is the main result of this theory that the growth of perturbations are damped by the Hubble expansion, which 
leads to a power law behavior for the time dependence of density perturbations. For example in the Einstein-de Sitter 
model which is believed to describe our Universe after recombination, perturbation amplitudes grow like (1 + z) . 
Only during the nonlinear stage with large density contrast does the evolution become faster. At the nonlinear stage, 
however, perturbations grow much faster, leading to the formation of gravitationally bound objects. 

The theory of linear density perturbations in a homogeneous medium was first developed by Jeans |90j|-|9l|. His 
study was motivated by the intention to explain the mechanism of star formation. We describe this theory below. 
First, however, its range of validity, i.e. the evolution of cosmological horizon, is discussed. 

The linear perturbations in the expanding homogeneous and isotropic Friedman Universe were studied by Lifshitz 
using a relativistic treatment. Relativistic theory, however, is necessary when the scale of perturbations is greater 
than the horizon, or when relativistic regimes are attained in matter condensation. In the most interesting cases 
such as perturbations in dark matter well inside the horizon after the equivalence epoch (when the energy densities 
of radiation and other components are equal), it is sufficient to consider nonrelativistic theory based on Newtonian 
gravity. Bonnor (see also was the first to study evolution of spherically symmetric perturbations in Newtonian 
cosmology. 

The theor y of linear density perturbations in the Newtonian treatment is developed in detail in various textbooks, 



M(r) oc r 



(21) 



III. GRAVITATIONAL INSTABILITY 



see e.g. 




A. 



Horizon scale and mass evolution 



The Newtonian treatment is only applicable on scales smaller than the horizon scale Xh = cH , The associated 
mass scale, defined as the mass contained within a sphere of radius Ajj/2, where H is the Hubble parameter, is given 




Figure 2: Two fractal structures a and b of dimension d ~ 1.6 generated by the same algorithm as c,d and e and endowed 
with the same lower cutoff, are compared and contrasted, a has no upper cutoff, while b has an upper cutoff at a distance d c , 
shown in e. In a the density decreases with distance as p oc r _1B , while in b, for r < d c the distribution is self-similar and 
equivalent to the one in a but for r > d c the distribution becomes homogeneous and p stays constant with distance (from [z3]). 



by 

MH = i np {¥) ■ (22) 

Beyond this scale events are causally disconnected and thus any correlation breaks down outside the horizon. Thus 
structures cannot form on scales larger than Xh- Mh monotonically increases with time because the distance that 
light travels increases with time. There are several regimes, separated by the moment of equivalence in energy densities 
of radiation and nonrelativistic matter: 

m » - { % : < z ^ 

Today the horizon scale is approximately 3000Mpc, which corresponds to a mass scale M ~ 10 22 M Q for an il = 1 
Universe. At recombination the total mass inside the horizon was then approximately (l/z rec ) -3 / 2 ~ 10 17 M©, where 
Mq = 2 10 30 kg is a solar mass. 



B. Self-gravitating ideal fluid: linear theory 



1. Fluid equations and background solutions 

Consider a perfect fluid with density p and pressure p in Euclidean space with a Cartesian (" physical" ) coordinate 
system j-j 1 . The fluid has a velocity field vf, the gravitational potential $ is induced by the mass density p distribution. 
All these quantities are related through the continuity, Euler and Poisson equations respectively. For a nonrelativistic 
fluid, i.e. for a fluid with p « pc 2 , they read |o3| [95l| |97| |9S||: 

^+di(pvi) = 0, (24) 



dvi 1 
— +v i B i v i +- 



+ VjdjVi + -d iP + = 0, (25) 



<9 2 $ - AirGp = 0, (26) 

where d 2 = didi. A cosmologically important solution of equations 124126( 1 is the one describing a spatially uniform 
fluid with zero pressure ('dust') on an expanding background [95j (99]: 

v?=H(t)n, (27) 



^+SH Po = 0, (28) 
Po=Po(t), (29) 

$o = l^Gpor 2 , (30) 
dH 2 4 

- + H* = - r Gpo, (31) 

where r 2 = r^i and all quantities depend only on time. 
In a comoving coordinate system, with 

r a = a(t)x a , (32) 
where a(t) is the scale factor, the relation between coordinate differences Ar^ and Ax a is 

dAr n dAx n da . dAx n , . . , 

~dT = a ~d~r + n - + H ®*°> ^ 

where 

H{t) = (34) 
a at 

Correspondingly, for velocity fields we have: 

v a (rp,t) = u a (xp,t) + Hr a = u a (xp,t) + u°. (35) 
Thus the solution ((27131(1 represents a uniform distribution of the fluid with zero peculiar velocity = and zero 



pressure po = 0. Pressure and density are linked through the equation of state p = p(p). The three equations I(24l2t 
together with the equation of state are a complete set, allowing one to study the temporal evolution of the density 
and velocity distributions as well as of the pressure and gravitational potential. 



Greek indices denote comoving coordinates, Latin indices denote physical coordinates, both take the values 1,2,3; the Einstein summation 
rule is adopted. 



2. Perturbed quantities 



As well known, solutions I|27I31|I represent an isotropic and homogeneous distribution of matter. In order to study 
density perturbations in the linear approximation assume that 

p(r i ,t)=p (t)[l + 5(r i ,t)], (36) 
v i (r h t)=v?(r j ,t)+Sv i (r j ,t), (37) 
$(ri,t) =# (ri,t) + *$(ri,t), (38) 



p(r l ,t) = Sp( n ,t), (39) 

where 6 = p ~^° . Here all perturbed quantities 5, Svi, dp and 5$ are assumed to be much smaller than the background 
quantities. All zero order values are given by (|27I31() . Assume also, that spatial together with temporal derivatives 
of perturbed quantities have the same order of magnitude as the quantities themselves. 

Note that it is not necessary for the smallness condition o n the perturbed quantities to hold over all of space. In 
particular, there could be a region in space where \Svi\ > \vf\ .100]. In this case the standard linearization procedure 
leads to different perturbation equations and consequently to different solutions representing a time dependent density 
contrast 5(ri, t). 

3. Linearized perturbations equations 

Rewriting 124I26[1 in comoving coordinates: 

d 1 1 

— + ZHp + -pd a u a + -u a d a p = 0, (40) 
at a a 

d 2 a du a 1 1 1 



.,, , Hu a + -updpiia H d a p+-d a $ = 0, (41) 

at at a ap a 

d 2 <P - AnGa 2 p = 0. (42) 

Here all quantities, except for H, d epend on the comoving coordinates 2 x a and the time t. Eqs. H40I42J) written in 
physical coordinates can be found in |100| for example. One arrives at the above results from I|24I26|I by using the 
transformation laws (d/dt) ph y S = {d/dt) com - Hx a d a and (d a ) phys = (l/a)(d a ) com . 

In order to obtain equations for the density contrast S in the linear approximation substitute I|36I39J) into equations 
(|40I42|) . Taking into account that spatial as well as temporal derivatives of perturbed quantities have the same order 
of smallness as the perturbed quantities themselves and using I|27I31|I , the perturbation equations read 

% + -d a 6u a = 0, (43) 
at a 

— ^ + H5u a + -d a S P + -d a 5<P = 0, (44) 
at a a 



2 If one also assumes that Hubble parameter can be disturbed (have spatial dependence) then the system of equations becomes overdefined. 
There is another approach llOll . however, where d a & and d a H are taken as independent variables in order to study density perturbations. 



d 2 c5$ - 4 7 rGa 2 / 9o<5 = 0, (45) 

where 6u a (xp,t) is first order quantity, because the unperturbed value is u° a (xp,t) = 0. 

The simplest way to find the equation governing density perturbations is to take the time derivative of Eq. I|43|) and 
use the divergence of Eq. (|44fl together with Eq. I|45|l . After rather long calculations one finds the final expression: 



d 2 S 85 - 

OP + dt a 



2H— - -§<9 2 <5 - 4irG Po 6 = 0, (46) 



where the relation 



is assumed. 



vl = (47) 
dp 



4- The Jeans criterion 

Eq. 1)46(1 governs the dynamics of density perturbations. It is a wave-like second order partial differential equation. 
Thus it is natural to introduce the Fourier transform 



5 = Y^h{t)e lk - x " (48) 



in order to split perturbations on different scales. 

Eq. H46J) can be rewritten in fc-space, taking into account that d a S — > ik a h: 



d 2 h nrr dh „ v 2 k 2 ., , , 

— = -2H- + (4.GP0 - ^)h, (49) 



where k a is the comoving wavevector and k — \/k a k a is the corresponding wavenumber. The comoving wavelength 
of the perturbative mode is given by / = 27r/fc, while the proper (physical) wavelength is simply A = al. 
The Jeans criterion (|49|l is governed by the wavelength 



where Aj separates gravitationally stable scales from unstable ones. Fluctuations on scales well above Aj grow via 
gravitational instability, while on scales smaller than Aj the pressure overwhelms gravity and perturbations do not 
grow. 

The first term on the right-hand side of l|49(l comes from the general expansion. In the static world initially 
considered by Jeans, such a term is absent, leading to the exponential growth of perturbations. In expanding space 
perturbations grow with time according to a power law. 

A very important quantity usually associated with the Jeans length is the Jeans mass (|50|l 

A^ = Ht) 3 ' (51) 
defined as the mass contained within a sphere of radius Aj/2, where p is density of the perturbed component. 



5. Multi-component system 

Perturbations for a given mode in a single component evolve according to (|49|l . When several components such as 
Cold Dark Matter (CDM), Hot Dark Matter (HDM), baryons and radiation are present simultaneously, it is possible 
to generalize (|46|l . Assuming gravitational interaction between components only, we arrive at 

d 2 hi nrr dh l „ v-^ , (v 2 )ik 2 , , , , 

~w = ~ 2H ii + {inGpo £ e ^ ~ ^h hi) > (52) 

3 

where the index i refers to the component under consideration, the sum is over all components and = Pi/^jpj- 
Notice that any smoothly distributed component (like the cosmological constant) does not contibute to the right-hand 
side of 



C. Applications 



Some important cases of matter content for the Unverse will be considered below. First we discuss perturbation 
dynamics in the dominant nonrelativistic component (baryonic or not). Second example is dark matter perturbations 
in the presence of a dominant radiation component. 

1. Einstein-de Sitter Universe 

First consider the dust dominated f2 = 1 Universe. This condition (17 is the density parameter) corresponds to a 
flat-type cosmological model, namely the Friedman solution of Einstein's General Relativity equations with curvature 
parameter k = 0. This model is thought to provide a good description of our Universe after recombination. To zero 
order a ~ t 2 / 3 , H = 2/3t and po — 1/QitGt 2 . For perturbations well inside the horizon we have 



,d 2 h 4 dh 2 
~d¥ + 3 ~dt ~ 3 



h = 0. (53) 



For modes well inside the horizon and still larger than the Jeans length the solution is 

h(k,t) = h i (k)(£\ 2/ \h 2 (k)(j- y \ . (54) 

As expected there are two solutions, one growing and one decaying. At late time, however, only the growing mode 
is important. Perturbations evolve proportionally to the scale factor or as (1 + z)~ l , where z is the redshift defined 
by: 

a{t) 

and ciq denotes the value of a scale factor today. 

Perturbations on scales smaller than the Jeans length cease to grow and oscillate with time. 

2. Mixture of radiation and dark matter 

Consider the radiation dominated Universe where a ~ t 1 ^ 2 and H = l/2t. The second component to be considered 
is a collisionless dark matter with vdm = 0. We still can use the Newtonian treatment on scales much smaller than 
the horizon size. 

Since the small scale photon distribution is smooth and the energy density is dominated by the radiation, the 
equation governing dark matter instability reduces to 

drhnM . dh DM _ ,.„.. 

+ = °- (56) 

This has the solution 

h DM (k, t) = h x (k) log (J^ + h 2 (k). (57) 

Perturbations in dark matt er component inside the horizon experience a slow logarithmic growth. This is the well 
known Meszaros effect |l02j . 

D. Initial spectrum of perturbations 

In the first example we have shown that perturbations on scales between horizon size and Jeans length A j <C A <C Xh 
grow asi5oc(l + z)~ 1 . In order to study perturbation dynamics one also needs to know the initial values of the 
perturbations at some moment in early Universe. 



One great possibility is provided by the CBR anisotropy measurements, because density inhomogeneities when 
photons were coupled to baryons can be extracted from temperature fluctuations observed in the CBR. Since 4?r — 
10 -5 it is usually assumed that at the moment of recombination S ~ 10 -4 [l5| . 

Perturbation amplitudes on different scales are usually represented by a power spectrum P(k), which is the Fourier 
transform of a previously introduced correlation function 96] 

«r) = ^3 / P(k)^4nk*dk. (58) 

There is no evidence that the initial spectrum contained any preferred scale, so it should be a featureless power law 

P(k) cx k n , (59) 

where the index n governs the balance between perturbation amplitudes on large and small scales. The value n = 
corresponds to white noise that has the same amplitude for every mass scale. The value n = 1 corresponds to a 
so-called Harrison-Zel'dovich scale invariant spectrum. Term 'scale invariance' means that perturbations had the 
same amplitude at the moment of horizon crossing. 



E. Damping of perturbations 

In addition to the Jeans scale some other cosmologically important scales appear in the theory of structure formation. 
They are related to physical processes that cannot be described within the perfect fluid approximation. However, 
fortunately such processes take place on limited scale intervals and outside such intervals the fluid description is still 
possible. We will discuss some dissipative effects, such as collisional damping of baryonic perturbations and free 
streaming of collisionless light particles. 



1. Silk damping 



Close to recombination the coupling between photons and baryons makes it possible for the former to erase pertur- 
bations of the latter. This is because at that time the free mean path of photons becomes larger, so they can travel 
from overdense into underdense regions dra ggin g baryons with them, thus smoothing inhomogeneities in the primeval 
plasma. This effect was discovered by Silk |103| . The physical scale associated with it is [TBI 

l s ~ 3.5(fi/i 2 )- 3 / 4 M P c, (60) 

which gives a mass scale 

M s ~ 6.2 10 12 (fWi 2 )~ 5/4 M Q . (61) 

This scale is close to the mass of a typical galaxy 1Q 1 Mq. However, Silk damping only affects baryonic perturbations. 
Moreover, it is important only around recombination when the coupling is still sufficiently strong to make photons 
drag baryons with them. 



2. Free streaming 

Another dissipative process is Landau damping or free streaming that originates from the free motion of collisionless 
particles on small scales. They can travel far if the velocity dispersion is large. This is important after particles 
decouple from the plasma and until they become nonrelativistic. 

The discovery of free stre aming was a dramatic moment for structure formation scenarios based upon Hot Dark 
Matter (HDM) models |l04j . The name Hot Dark Matter means that particles composing such matter were ultrarel- 
ativistic at equivalence. Thus their velocity dispersion was near the speed of light c. The maximum distance scale 
travelled by collisionless particles from decoupling can be estimated to be j9?| 

l FS * 0.5(^r 4 / 3 (fW 2 ) 1/3 Mpc, (62) 



and a corresponding mass scale is of the order of superclusters of galaxies or even larger if the particle mass is 
itldm < 30 eV. 



Cold thermal relics which compose the Cold Dark Matter (CDM) are slow enough so that free streaming can be 
neglected on cosmologically important scales. Therefore Landau damping affects only light particles like neutrinos 
with m v ~ 10 eV, and in the Universe dominated by HDM, all perturbations on scales smaller than superclusters of 
galaxies are erased. At the same time, after particles become nonrelativistic at z nr , their velocity dispersion becomes 
small enough to make free streaming negligible. 

F. Structure formation at late times 

1. Nonlinear clustering 

In previous sections we dealt with the linear evolution of cosmological perturbations. From long after recombination 
and even earlier until almost recent times, such a treatment is justified to describe the growth of inhomogeneity in 
the Universe because the condition 8 <C 1 is sutisfied. In recent times, say at z ~ 10, the nonlinear behaviour of 
perturbations becomes important. During the nonlinear stage gravitationally bound objects such as galaxies form. 
Nonlinear evolution is a rapid process, where not only gravitational effects are important. In particular during the 
formation of galaxies, various dissipation and relaxation processes take place [9(|. 

We will not go into the details of galaxy formation here. The LSS formation is the subject of this section. Here 
the main interaction remains gravity. However, the theoretical description based on linear equations for an ideal fluid 
becomes inadequate. 

Usually N-body simulations are employed to study nonlinear clustering |96(. However, some useful simplified models 
are still possible even in the nonlinear stage, because numerical simulations provide limited physical insight into the 
physics of gr avitational clustering. Among nonlinear approximations, the most famous are the Zel'dovich approxi- 
mation |l05j and the spherical collapse model (see for example |&|). The key point in the Zel'dovich model is that 
during collapse in an almost spherically symmetric overdense region, gravitational interaction amplifies asymmetry. 
Therefore, the final structure acquires a preferred direction and the final collapsed body will look like a 'pancake'. In 
the spherical collapse model, on the other hand, it is assumed that spherical symmetry remains valid during the entire 
period of collapse. It allows the splitting of the spherical overdense region into concentric shells and the evolution of 
each shell can be studied separately, which sufficiently simplifies the problem. This model will be described in detail 
in the next chapter. 

2. Structure formation scenarios 

Historicall y tw o different pictures of structure formation were considered, namely the HDM (see |10(t| ^ and CDM 
models (see |lQ7j V We will discuss them briefly below. 

3. HDM models 

The neutrino dominated Universe with m v ~ 10 eV is a typical HDM model. The HDM model is associated with 
so-called "top-down" scenario, where structures form on large scales first. This is so because the Jeans mass for HDM 
is of the order of the supercluster mass or even higher. At the same time, free streaming erases perturbations on 
smaller scales. Thus only when perturbations reach the nonlinear regime on large scales can they induce fragmentation 
on smaller scales. 

Usually it is assumed that large scale perturbations become nonspherical according to the Zel'dovich model and 
thus the LSS look like a "net" of density condensations separated by huge voids. Simulations agree with such a 
picture. 

The HDM model is in good agreement with observational data on scales larger than lOMpc. On smaller scales, 
however, HDM simulations can agree with observed the correlation function of galaxies only if the epoch of pancaking 
takes place at z ~ 1 or less, which is too late, because we can see galaxies and quasars with much greater z. 

The crucial cosmological property of HDM with neutrinos, as was mentioned above, is the damping of perturbations 
on small scales due to free streaming. Neutrino dominated models (Sl„ ~ 1) alone could not describe the real Universe 
because on scales smaller than ~ 100 Mpc no structure appears at all. 

One important prediction of HDM models with neutrinos is the existence of large smooth halos around galaxies. At 
the end of collapse, during formation of galaxies, the baryonic component can dissipate its energy via collisions, but 
the neutrino component cannot. Thus neutrinos remain less condensed than baryons, forming large galactic halos. 



4. CDM models 



CDM models do not have trouble with free streaming, because their particles have negligible velocities at decoupling. 
Moreover, the Jeans mass for typical CDM model lays well below 10 6 M Q . Thus, perturbations start to develop on 
small scales simultaneously with perturbations on large scales. 

In the CDM models the important feature is the weak growth experienced by perturbations between horizon crossing 
and equivalence (see flIII C 2f) . This means that the density contrast increases when we move to smaller scales, or that 
the perturbations spectrum has more small-scale power. 

After collapse the first structures in CDM models virialize through violent relaxation |108| into gravitationally 
bound objects that form galactic halos. Structures form in a self-similar manner from small to large scales, in other 
words according to a 'bottom-up' scenario. 

Pure CDM models, however, fail to predict the observed correlation function for galaxies on large scales. If one 
wants to retain the CDM hypothesis, the simplest way is to reduce the matter density. This shifts matter-radiation 
equivalence to a later epoch, resulting in redistibution of power in the spectrum of perturbations in favour of larger 
scales. 

Today the ACDM model with f2 tot = 1 and = 0.7 is considered to be the best fit to the full set of observational 
data (see E.W. Kolb, in these proceedings Q). 



IV. NEUTRINOS AND STRUCTURE FORMATION 



In previous chapter we described the evolution of perturbations and we saw that the nature of dark matter particles 
is crucial to determining the way structure formation develops. In spite of the fact that a lot of candidates for CDM 
particles exist, there is presently no experimental evidence of such particles. On the other hand, neutrinos are the 
only candidates for DM known to e xist. 

'Light' neutrinos (m„ <C lMeV) 109], namely neutrinos that decouple while still in their ultrarelativistic regime 
(see below), may provide a significant contribution to the energy density of the Universe ( fl„ ~ 1). Models with light 
neutrinos were extensively studied in the eighties; a large literature exists on this subject [llOj . 

The key prediction of the cosmological model with neutrinos is a cellular structure on large scales (see Fig. 0. 
The qualitative drawing of cellular structure of the Universe is represented in Fig. Ruffini and collaborators have 
studied such models with particular attention to the problem of clustering on large scales and its relation to the fractal 
distribution of matter. In the following, we will outline some of these ideas. 



A. Neutrino decoupling 

The cosmological evolution of a gas of particles can be split into two very different regimes. At early times, the 
particles are in thermal equilibrium with the cosmological plasma; this corresponds to the situation in which the rate 
r =< avn > of the reactions supposed to mantain the equilibrium (such as v e + v e «-> e + + e~ <-> 27 in the case of 
electronic neutrinos) is much greater than the expansion rate, given by the Hubble parameter. The gas then evolves 
through a sequence of thermodynamic equilibrium states, described by the usual Fermi-Dirac statistics: 

f(P) = ut?( \ 1 wi tTTT ' ( 63 ) 

exp [(E(p) - fj.)/k B T\ + 1 

where p, /1 and T are the momentum, chemical potential and temperature of neutrinos respectively, and ks is the 
Boltzmann constant. 

However, as the Universe expands and cools, the collision rate Y becomes smaller than the expansion rate; this 
means that the mean free path is greater than the Hubble radius, so we can consider the gas to be expanding without 
collisions. It is customary to describe the transition beetween the two regimes by saying that the gas has decoupled 
from the cosmological plasma. 



1. The redshifted statistics 



Since in a spatially homogeneous and isotropic Universe described by the Robertson- Walker metric, the product of 
the three-momentum p(t) of a free particle times the scale factor a(t) is a constant of the motion: 

p(t) ■ a(t) = const, (64) 




Figure 3: Qualitative drawing of the evolution of the cellular structure of the Universe. The vertical axis represents decreasing 
values of the redshift. z — represents the cosmological observer today, embedded in his own 'elementary cell'; z\ z%, z n 
are the redshifts of the centers of the successive 'elementary cells' seen by this observer, z = z 1 represents the surface of last 
scattering (decoupling) of the CBR. 6 is the angular size subtended today by an 'elementary cell' at the time of decoupling. 
Finally, z = z nr is the time at which the cells formed via gravitational instability. 



each particle in the gas changes its momentum according to this relation. This fact, together with Liouville's theorem, 
implies that the distribuition function after the decoupling time t d (defined as the time at which T — H) is given by 



f( P ,t>t d ) = f 



aft) 
a d 



P,td) = 



1 



exp 



(65) 



where the subscript d denotes quantities evaluated at the decoupling time. 

Now we turn our attention to the special case of neutrinos with m v < 10 eV. The ratio T/H, as a function of the 
cosmological temperature, can be evaluated using quantum field theory [l5j 



r 

H 



T 



lMeV 



(66) 



as long as T > m. Therefore, neutrinos decouple from the cosmological plasma when T = T d — 1 MeV. Since 
kT d 3> tog 2 , many of the particles satisfy pc 3> mc 2 and then, when performing the integration over the distribution 
function l|65(l . we can safely approximate: 



f( P> t>t d ) = f(^p, t d 



1 



exp 



a d 



pc- fid) /k B T d 



(67) 



+ 1 



since the tail of the distribution function for which mc 2 3> pc gives little contribution. 

In the following, we compute the mean value of physical quantities over this distribution. It will be useful to 
consider two limiting regimes, namely the nonrelativistic one and the ultrarelativistic one. They correspond to two 



approximations for the single particle energy |7! 



E~mc 2 kT<^mc 2 , NR, 
E~pc kT^>mc 2 , UR. 



(68) 



We emphasize the fact that this substitution must be performed only in the function to be integrated, and not in the 
distribution function. The approximation l|67|) depends only on the fact that the particles are ultrarelativistic at the 
time of decoupling, and then it is valid even when kT <C mc 2 . 

Then, with a suitable substitution of variables, all the relevant integrals can be recast into a very simple dimen- 
sionless form: 

r n (0= [°° — .. ynd * - - , (69) 

Jo ex P i(y ~ 0} + 1 

where £ = [i^jkT^ is the dimensionless che mica l potential or degeneracy parameter. These integrals can be expressed 
using Riemann zeta and related functions 

2. Energy density of neutrinos 

The present density parameter of neutrinos can be easily evaluated using the method outlined in the previous 
section. The energy density is given by: 

/>OC 

P»+»(to) = -^r E(p)f(p,t )d 3 p, (70) 



density p c = 1.054 h 2 ■ 10 4 -^4, we obtain 



l p Jo 

where g is the number of helicity states and hp is Plane's constant. By normalization with respect to the critical 

Sl^h 2 ^ 1.10 -lQ^g^A^), (71) 



and hp : 
[H |Z3: 



where A(£) is defined as follows 

h(0 + h(-0 



2/ 2 (0) 477(3) 



-\e+ Mm\ +4^(-i) fc+i e ^ 



k 3 

k=l 



(72) 



and rj(n) is the Riemann eta function of index n. 

The term l2{— £) appears because we have to take into consideration the presence of antiparticles, for which the 
relation & = — £ v holds. This result follows from the fact that if we consider a reaction such as 

v + u< — ► ... < — >7 + 7, (73) 

then since the chemical potentials of the initial and final states have to be equal and the chemical potential of the 
latter is equal to zero, it follows that £ p = — £„. 



3. Recent constraints on the neutrino mass m„ and degeneracy parameter £„ 

We discuss below very briefly recent bounds on the chemical potential and mass of neutrinos |l!2j that can be used 
to compute the recent bound on O^+p. 

a. Neutrino mass A recent laboratory limit on the electron neutrino mass comes from tritium j3 decay |22l |. 
These data give limits 

m„ e <2.5eV. (74) 

At the same time, no direct measurements or constraints on muonic and tauonic neutrino masses exist. Moreover, 
it is still unknown whether neutrinos are Majorana or Dirac particles. Very recent data from neutrinoless double 
decay |l!3j give also lower bound on Majorana mass: 



(0.05 < m„ cc < 0.86) eV. 



(75) 



b. Chemical potential The first constraints on the neutrino degeneracy parameter from cosmological nucleosyn- 
thesis were obtained in 

El- 

It was shown later |l6j ] that a small value of £ e coupled with large values of |£/*,t| can 
lead to cosmological nucleosynthesis abundances which are consistent with observations. It is found in particular that 

< £e < 1.5, (76) 

with the additional constraint F(^) + F(£ T ) ps F(10£ e ), where F(£) = £ 2 + £ 4 /27r 2 . In particular this implies 

l^,r|<10£e. 

Recent data both fr om c osmological nucleosynthesis and CMBR |l!8| strongly constrain neutrino degeneracy pa- 
rameters. Orito et al. |l!5| give surprisingly wide const r aints , £ e < 1.4 and \£n, T \ < 40. Other papers give essentially 
stronger constraints using additional assumptions |l!6j |- [TT^ | 

£e < 0.3 

l^.rl < 2.6. V 7 > 

c. Neutrino oscillations When one consider different chemical potentials for all n eutri no flavors at the epoch 
prior to cosmological nucl eosyn thesis, neutrino oscillations equalize chemical potentials |119| if there is enough time 
for the relaxation process 120]. On the basis of large mixing angle solution of the solar neutrino problem which is 
favored by recen t dat a |1 2 1| , cosmological nucleosynthesis considerations constrain the degeneracy parameters of all 
neutrino flavors Il2 



|£| < 0.07. (78) 

However, the situation when flavor equilibrium is not achieved before cosmological nucleosynthesis is also possible. 
Thus in the following we consider quite high values of the degeneracy parameter and assume its value positive without 
loss of generality. 

The main result that comes from considering oscillations is that masses of different neutrino species are nearly 
equal: m Ve ~ m v ~ m„ T . 



One can see that quite high values of the neutrino energy density are still possible if we assume the recent constraints 
discussed above. In particular, if one consider two Dir ac ne utrino flavors with equal masses and chemical potentials 
(y e gives very small contribution to fi„), £ M = £ r < 2.6 |ll7j | and m VlA — m Ur < 2.5 eV j^, one gets the upper bound 

fl v+j? < 0.45. (79) 



4- The Jeans mass of neutrinos 



In neutrino dominated Universe the first possible structure occurs when these particles become nonrelativistic, since 
at earlier times free streaming erases all perturbations. At this epoch the cosmological redshift has the value lljj 



1.698 10' 



where 



m = — W) — 



48?7(5) 



;f + 8r)(2)e + 4877(4)^ + 48^(-l)" +1 - 



n=l 



(80) 



(81) 



The basic mechanism of fragmentation of the initial inhomogcncitics in an expanding Universe is the Jeans instability 
described in the previous section. However, in the calculation of the Jeans length of nonrelativistic collisionlcss 
neutrinos, we cannot use the velocity of sound obtained from the classical formula (|47|l . In fact, since the particles are 
collisionless, their effective pressure is zero and this would lead to a vanishing Jeans length, meaning that even the 
smallest perturbation would be unstable. This is not the case since in the absence of pressure, another mechanism 
works against gravitational collapse, namely the free streaming of particles (see 9IIIE2JI . The characteristic velocity 
associated with this process is simply the dispersion velocity y/ < v 2 > /3, where the factor 3 comes from averaging 



over spatial directions. Thus, we have to make the substitution 
< v 2 > can be obtained using the method described above: 



>< v > /3 |2J|. The correct expression for 



2 *\ 



(82) 



where T v q = 1.97 K is the present temperature of neutrinos. 

As a result, the Jeans mass grows in the UR regime and decreases in the NR regime [l2j|. The evolution of the 
Jeans mass of neutrinos for m„ = 2.5 eV and £ = 2.5 with redshift z is described by Fig. 0] It is clear that for such 
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Figure 4: The Jeans mass dependence on redshift for neutrinos with mass m„ = 2.5eV and degeneracy parameter £ = 2.5. 



values of the neutrino mass the peak of the Jeans mass lies above 10 Mq and the corresponding comoving Jeans 
length is Ao > 100 Mpc. On the other hand, the value of the Jeans mass today is still larger than the mass of massive 
galaxy 10 12 M Q . 

Finally, the maximum value of Jeans mass at the moment (|80|l is 24] 



Mj(z nr ) = 1.475 10 17 M q5 „ *N V *(—!!-)-* A{£)~* 

lOeF 



(83) 



The peak of the Jeans mass as a function of the degeneracy parameter for different fixed values of the energy density 
as well as with constant mass ra v = 2.5 eV is shown at Fig0 

By comparing different curves with a fixed value of £ one can find the well known result that the Jeans mass 
increases with decreasing of neutrino mass. With the growth of degeneracy parameter, however, the neutrino mass 
decreases in the beginning, and its different values correspond to different points at the same curve. 

The space above the dashed line at Fig. represents the region in which the neutrino mass is less than 2.5 eV. It 
iss interesting to note that this value of m u is still sufficient to get Sly = 1 with £ sa 4. 



B. Subsequent fragmentation model 



In this section we will describe a model of structure formation which explains the observed fractal distribution 
matter. The key point of this model is the existence of upper and lower cutoffs in the fractal. The upper cutoff 
appears due to causality of the elementary cell, and the lower cutoff corresponds to the time when the dark matter 
ceases to dominate in the formation of the structures 621. 
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Figure 5: The Jeans mass dependence on the degeneracy parameter with a fixed value of energy density, curves (1-4). Curve 
(1) corresponds to energy density Q. v — 0.11. Curve (2) corresponds to fl v = 0.3. Curve (3) represents the neutrino energy 
density Q. v = 0.5 and finally curve (4) gives Jeans mass for Q. v — 1. The dashed line represents Jeans mass dependence on the 
degeneracy parameter with fixed neutrino mass m„ = 2.5 eV. 



1. Nonlinear model of spherical collapse 

Following Ruffini et al. : 8l] we consider the first spherical perturbation in the Friedmann Universe. For simplicity 
we neglect the interaction between neighboring perturbations, and treat a given elementary cell formed at the epoch 
z nr as a spherically symmetric region. If the mass of an elementary cell is sufficiently larger than the Jeans mass, 
the mass density will dominate the pressure p <C p and then the sphere will expand freely from the expansion of the 
Friedmann background Universe. Then the equations governing the dynamical evolution of this shell can be obtained 
by solving the Einstein field equations in the Friedmann Universe in the metric 

ds 2 = -dt 2 + e 2A(x,t) dx 2 + r 2 {x,t){d0 2 + sin 2 0# 2 ), (84) 

where % is a comoving radial coordinate and r(x, t) is the radius of a 2-sphere. This metric is very similar to the 
Robertson- Walker one, because the spherical overdense (underdense) region behaves like a closed (open) sub-universe. 
The solution of the Einstein equations with the above metric gives an equation for the total energy of a test particle: 

. 2 2Gmix , . , . 

n hei(x) = 0, (85) 

n 

where mi(x) an d £i(x) are respectively the gravitational mass and the negative total energy of the test particle at 
radius r% of a 2-sphere representing the first perturbed region. 

We assume that £i(x) is constant in time since there is no dissipation and that at the moment of separation from 
the expansion of the background the expansion rates are the same (a/ a )i = ( 7 V r )i- If l ne mean mass density of the 



sphere is defined by 



Pi(t) 



_ mi(x) 
" fr?(x,t)' 



(86) 



then we can write the energy ei(x) as 



ei(x) = yGM*i)-R?*i 



(87) 



Here we have assumed a Friedmann Universe with k = as a background. i?i = r(x,ii), pi(ii), and pu{ti) are 
respectively the radius of the sphere, its mean density, and the mean density of the Universe at the epoch of separation 
t = t\. Furthermore, 8\ = p\{t\) / ' pu{ti) — 1 and H\ is the Hubble parameter at z = z\. Eq. (|85|) can be solved in 
parametric form: r\ = ri(Oi), t — t(0i), where the parameter 0i is a conformal time which has an initial value 9\ 
related to the initial value of the density contrast 5\ by 



1 



1 + Si - 2eSi 
l + 6i 



where e = —1, +1 corresponds to underdense and overdense regions respectively, being the analog of the curvature 
parameter k in the usual Robertson- Walker metric. Using the relation between cosmic time t and redshift z in a flat 
(k = 0) Friedmann Universe, we arrive at the following expression for z as a function of 0: 



1 + z = (1 + zx) 



l + 6x 



3/2 



(9i - sinGi - 6»i + sinfli) 



-2/3 



(89) 



At the same time, an expression for <5i can be obtained: 

3 



*i(t) = 8 
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l + z 
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(l + ^i) ; 



-(l-cos6i)- 



(90) 



This set of equations determines the evolution of the elementary cell. 



2. Successive fragmentation 

After the separation of the elementary cell from the Hubble flow, the Jeans mass will continue to fall and perturba- 
tions of smaller masses will be able to detach themselves from the expansion either of the cosmological background or 
of the larger parent elementary cell. The perturbed regions in the background will follow the evolution given above, 
but the perturbed regions in the parent elementary cell clearly follow a different evolution. 

Here we discuss the simplified case of a spherically symmetric pertubation separating from the expansion flow of 
its parent cell. We can proceed exactly as we did for the initial elementary cell perturbations, remembering, however, 
that now the perturbation detaches itself from the still expanding parent cell and not from the Hubble flow. Thus 
the density and evolutionary state of the parent cell, and not the one of the Friedmann background, will enter into 
our calculation of the evolution of successive perturbations. 

For the perturbation of n we thus have 

s »(^< i TTf) 3 Trw ,1 - C08e '' r3 - 1 ' (91) 

where A„ is a parameter in which the whole history of the previous fragmentation is summarized. It has the expression 

A -1 + (4-l+An-l) 3 Si^Q^-l rq9 s 

" \i + z n ) a+s^r (l-cose^)*' 

where Q^—i ^ s the value of conformal time when the n-th fragmentation appears inside the (n— l)-th one, and Ai = 
for the first fragmentation. 
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Figure 6: Qualitative illustration of subsequent fragmentation process in the spherical model. The final picture obtained from 
such a mechanism looks like a typical fractal. 



3. The fractal model 

To explain in a simple way the mechanism of this model, we start with a simple twofold scenario in which each 
condensation gives birth to two daughter condensations as soon as the value of the Jeans mass of the parent cell drops 
to half of the initial value. In short, we take a condensation of mass Mj, then reach the redshift at which the Jeans 
mass has become Mj/2. At that epoch inside the initial condensation two new daughter condensations originate, 
each of mass Mj/2. As we continue this process, we will have four granddaughter condensations, then eight, sixteen 
and so on. This process can easily be generalized to the occurence of N fragments at each step. 

Since the mass of the initial condensation as well as the behaviour of the Jeans mass and length are given, the 
only free parameter at each successive step is the amplitude of the perturbation S n . Our goal is to reproduce at the 
end of the fragmentation process the expected fractal distribution and, therefore, to select at each step the suitable 
perturbation 5 n for this purpose. It has been shown [jj that such a process leads to a system with fractal dimension 
Dp given by 

Ioe N 

Df = T^y, (93) 
log A 

where A = rj_i/rj is a constant. In our case we assume Dp = 1.2 H3. 

In the spherical model dccribed above the density need not be uniformly distributed inside the radius r n : any 
spherically symmetric perturbation will clearly evolve at a given radius r n in the same way as a uniform sphere 
containing the same mass. We assume then without loss of generality that 

= HZ (94) 

The process of fragmentation is shown qualitatively in Fig. |H| This picture is reminiscent of a fractal. Note that 
the position of each new fragment inside previous one is not important. Moreover, predictions of the model are quite 
insensitive to the number N. The key point is that even for a random number of fragments appearing at each step 
the resulting density distribution is still a fractal. One difficulty with this model is that it predicts too large a density 
contrast So today at the galactic scale. 

In order to avoid such high values of the density contrast the authors introduced a suitable lagging time factor r. 
This factor is a function of N that introduces a time delay in the formation of each daughter condensation. It was 
shown that a valid phenomenological relation leading to realistic values of density contrast today is 

N 2 3 

t(N) = — + ^N. (95) 

Clearly, the existence of the lagging factor r is not in contradiction with the Jeans instability picture and only means 
that the fragmentation occurs somewhat later than the time at which the necessary condition is fulfilled. 

In the model under consideration the only free parameters are the initial density contrasts S n at every step. These 
are chosen in such a way that Eq. (|93|l is satisfied at each step. Thus the natural question arises, what is the form of 
the initial spectrum? It can be obtained by following backward in time the evolution of perturbations. The result is 



quite surprising and simple. The spectral index of initial spectrum at z nr (see <|59|) ) is n = 0, which corresponds to 
white noise. Note also that agreement with the observed CBR anisotropy can be obtained within the framework of 
this model [8fi| . 
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